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Abstract. Wc demonstrate that a stable particle with flavor mixing, confined 
in a gravitational potential can gradually and irreversibly escape - or "evaporate" 
- from it. This effect is due to mass cigenstate conversions which occur 
in interactions (scattering) of mass states with other particles even when the 
energy exchange between them is vanishing. The evaporation and conversion 
are quantum effects not related to flavor oscillations, particle decay, quantum 
tunneling or other well-known processes. Apart from their profound academic 
interest, these effects should have tremendous implications for cosmology, e.g., 
(1) the cosmic neutrino background distortion is predicted and (2) the softening 
of central cusps in dark matter halos and smearing out or destruction of dwarf 
halos were suggested. 



PACS numbers: 03.65.-w, 14.60.Pq, 14.80.-j, 95.35. +d 

1. Introduction 

What happens to a particle trapped in a gravitational potential? The answer is 
obvious: it will remain gravitationally bound forever. The particle will never escape 
from the potential even if it interacts with other particles, provided these interactions 
preserve the total energy of the particle. Paradoxically, this answer turns out to be 
false for a flavor-mixed particle, such as a neutrino, for example. In this paper we 
show that a stable mixed particle, which bounces in a gravitational potential and 
scatters off other particles from time to time, gradually escapes — or "evaporates" - 
from it. Strictly speaking, the probability to detect the particle inside the potential 
decreases monotonically with time and, of course, the probability of its detection 
elsewhere increases. Such 'evaporation' is a result of the conversion effect of mass 
eigenstates (in particular, of a heavier state into a lighter one) — another process 
that has not been addressed. We underscore that particle evaporation and mass- 
state conversion processes have nothing to do with flavor oscillations, particle decay 
or quantum tunneling. 

Besides its profound academic and general physics interest, the evaporation 
effect shall have tremendous implications for cosmology. Almost one quarter of 
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the energy and matter in the universe is dark matter, whose presence is revealed 
via its gravitational interaction only (see, e.g., Ref. pQ, for an overview), e.g., via 
gravitational strong and weak lensing, via stellar and gas kinematics, etc. The analysis 
of how the observed large-scale structure of the universe has formed supports the cold 
dark matter model with a cosmological constant (ACDM), in which non-relativistic 
(i.e., 'cold') elementary particles gravitationally collapse to form dark matter filaments 
and halos. Gravity of the halos traps normal gas, which either collapses further to 
form stars and galaxies, or remains hot in galaxy clusters where it is detected with 
X-ray observations. Usually, dark matter dominates the mass of a self-gravitating 
system, e.g., the gas mass in clusters is only about 10% of the total and the rest 90% 
is dark matter. The origin of dark matter remains unknown, but the currently favored 
models, e.g., the axion dark matter and weakly interacting massive particles (WIMPs), 
admit them to be mixed particles: an axion can be mixed with a photon and WIPMs 
are a mixture of several supersymmetric particles (see, e.g., Ref. [H 03 21 |S] , f° r recen t 
reviews) . How gravitational collapse of such flavor- mixed CDM occurs in the presence 
of evaporation and to what extent the formed dark halos and the entire large-scale 
structure of the universe are affected by this process is of paramount importance for 
modern precision cosmology. We suggest that dark matter evaporation may be able 
to address the outstanding problems of the CDM model: [H [7] the overproduction of 
small-scale halos known as the 'missing satellite problem' and the CDM prediction of 
diverging density profiles (cusps) in halo cores, not observed in dark-matter-dominated 
galaxies. 

After nearly eighty years since its discovery and decades of observations and active 
search, the nature and properties of dark matter remain unknown. The possibility to 
constrain its mixing properties is very exciting; the study can readily be done with 
cosmological simulations. 

Another intriguing cosmological implication deals with the cosmic neutrino 
background (CNB) — the copious relic neutrino leftover from the big bang epoch [8], 
whose present day temperature is < 2 K. Although neutrino masses are unknown, an 
estimate [H] of m ~ (Am 2 ) 1 / 2 ~ 0.1 — 0.01 eV implies that cosmological neutrinos are 
nonrelativistic and can be trapped in CDM halos with escape velocities of the order of 
a few thousand km s _1 , so they can be affected by conversions and evaporation. Thus, 
if one would be able to detect CNB on Earth, one can notice that its flavor composition 
is noticeably skewed compared to the traditionally expected equipartition distribution. 

2. Description 

The propagation (mass) and interaction (flavor) eigenstates of mixed particles are 
related by a unitary transformation, = J^j^ij l m j)> where |/) and \m) denote 
the flavor and mass states, and U is a unitary matrix. For simplicity, we consider 
a stable two-flavor particle, so U is a 2 x 2 rotation matrix through an angle 8 (8 
being the mixing angle). The particle interacts in the flavor basis, so the interaction 
matrix is diagonal in this basis, V = diag(V" Q , Vp) : but it is non-diagonal in the mass 
basis, V = WVU. The masses of the mass-states are mn and mi < m^, i.e., heavy 
and light. In general, mass states have different four-momenta and propagate along 
different geodesies. 

The process of evaporation is easily explained with a specific example of an 
electron- muon neutrino system in one dimension. Let's consider a thought experiment 
in which an electron neutrino is created inside a gravitational potential, 4>{x). Energies 
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and momenta of the mass states depend on the production process. By choosing them 
appropriately, we make the velocity of the light state to exceed the escape velocity of 
the potential, v esc , and the velocity of the heavy state to be smaller than u csc . Then, 
the trajectory of \rrih) is bound inside (j>(x) and that of \mi) is unbound and escapes to 
infinity (e.g., it is always so if |mj) is relativistic or massless). In time, only the heavy 
state is left inside and if nothing else happens, this is the final state of the system. 

Next, we let \rrih) to scatter off another particle (mixed or not) of mass M. 
Scattering of particles occurs through V, which is non-diagonal in the mass basis. 
The off-diagonal terms, V\h = Vhi = (V e — V^) cos 9 sin 9, couple different mass states, 
so that the amplitudes of \mh) and \m{) can be non-zero upon scattering, i.e., the 
scattered wave is a mixture of mass states again. Hence, there are two channels 
m/j — » rrih and mh — > mi, which are the elastic scattering of a mass-state and 
conversion of a larger-mass state into a lower-mass one. 

The total energy and momentum must be conserved in both processes. All is 
trivial for elastic scattering, so we consider the conversion process only. In the center 
of mass frame the momentum and energy conservations are Ph + Pm = = p\ + p' M 
and (ml + pl) 1/2 + (M 2 + P 2 M ) 1/2 = {m 2 + p' 2 ) 1 / 2 + (M 2 + p' M 2 ) 1 / 2 , where 'prime' 
means 'after scattering' and c = K = 1 throughout, that is 

{ml+pl) 1 ' 2 + {M 2 +p 2 h y' 2 = (mf+pf) 1 / 2 + (M 2 + P ' 2 y/ 2 . (1) 

Two special cases worth consideration: (a) M — > oo and (b) mh M. The first 

corresponds to the interaction with an external potential (an axion in a -B-field) or with 
a heavy particle (a neutrino interacting with matter) . The second is interesting for the 
self-interacting WIMP model (Am = m^ — mi <C mi ensures that both supersymmetric 
particles are stable, see below). In both cases, we have m 2 + p\k m 2 +p\ , i.e., the 
kinetic energy of the secondary exceeds that of the primary. (Note that the reverse 
conversion process, mi — ¥ mh can, therefore, be kinematically suppressed.) In case (a), 
\m{) can be relativistic with p[ ~ m^, if 3> mi, or non-relativistic with the velocity 
v[ ~ (2Am/m;) 1 / 2 , if Am -C mi; in case (b) both mass states are nonrelativistic and 
v[ - (2 Am /mi) 1 / 2 . 

After the scattering, the two mass states propagate along the geodesies again, and 
if v'i > v csc , then \m{) leaves the potential whereas \mi) remains trapped. Therefore, 
upon such a process the amplitude of the heavy state decreased irreversibly. The total 
probability is still unity, but the probability to detect a particle (an electron neutrino, 
for example) inside the potential has decreased and the probability of its detection 
somewhere outside has become larger. Of course, the overall energy is conserved: 
the light state climbs up the potential and looses energy (and a massless particle is 
redshifted). By repeating this cycle, one can further decrease the amplitude of the 
trapped state; colloquially speaking, the particle "evaporates" from the potential well. 

We illustrate the effect numerically. We set M — > oo (hence, the velocity of the 
scatterer is negligible) and Am -C m; (hence, all particles are non-relativistic). The 
conservation of four- momentum, Eq. ([!]), in this case is 

p 2 h /2m h = p 2 /2mi - Am. (2) 

To study the evolution of the non-relativistic mass states, we solve the two-component 
Schrodinger equation, 

id t \m(t, x)) = (H iicc + + V) \m(t, x)) , (3) 

where \m(t,x)) = (ah(t,x),ai(t,x)) is the state vector. Here the free particle 
Hamiltonian, H llcc — di&g(—d 2 x /2mh,—d 2 x /2mi — Am), satisfies condition (pi). 
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Figure 1. Space-time diagram of a flavor-mixed particle in a gravitational 
potential obtained by solving a two-component Schrodinger equation. The 
probability densities of the light (blue) and heavy (orange) mass eigenstates are 
shown. The potential is localized between x ~ —1 and +1, scattering occurs at 
x ~ 0.1. At t = a particle is created in a flavor state. The first scattering occurs 
at t ~ 3.5. Each scattering produces forward and reflected wave packets; those 
corresponding to light mass eigenstates escape and heavy remain bound. 

Gravity enters via H SIav = diag(m/j0(a;), mi<f>(x)), where we chose a model attractive 
potential with exponential screening, cf>(x) — (f>oe~( x ' x «> (l + (x/x g ) 2 )^ 1 , where 4>o < 
determines its depth and x g sets its size (x g ~ 0.5 in computational units). For 
simplicity, we set — and V e ^ 0. The scatterer is placed off-center, at 
x s ~ 0.1 and the scattering potential was chosen to be well- localized in space, 
V e = Voe-^-^V*.) 2 ^ + (( x - x s )/x v f)- x , where V < and x v ~ 0.005; the 
actual shape of V(x) does not affect qualitative results so long as x v is sufficiently 
small. We solve the initial value problem with gaussian wave packets and different 
momenta of the mass states at t = 0. 

The exact numerical solution of the Schrodinger equation of a mixed particle is 
shown in Fig. [T] It represents the space-time diagram of the probability density 
(i.e., the amplitude squared) of the heavy (orange) and light (blue) mass states. 
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(The amplitudes of flavor states can be obtained easily as a linear combination, 
\fi) — J^jUij Initially, the particle is in a flavor state. This is a coherent 

combination of two mass states, which propagate differently and, in time, the light 
state escapes. In contrast, the heavy state bounces off the potential and at t ~ 3.5 
scatters off the small-scale potential at x = 0.1 for the first time. At this time, 
forward-scattered and reflected waves emerge, both arc mixtures of two mass states. 
The light state escapes to infinity again but the heavy one keeps bouncing. Repetitive 
scatterings produce more outgoing light state wave packets at later times. 

The purpose of this toy model solution was to demonstrate the existence of the 
evaporation and conversion processes. We believe it is a very beautiful result on its 
own. Applications to real physical systems may require modifications. For instance, for 
dark matter particles, the gravitational potential shall be computed self-consistently 
from the particle distribution. 

3. Extension and remarks 

Gravitational field is efficient but not the only way to separate mass eigenstates. 
Other possibilities include scattering with gravitons (a rather weak, but not impossible 
process) or just ballistic propagation with different velocities. Thus, the above model 
of a gravitationally confined particle can be reduced to a simpler "particle-in-a-box" 
system, whose walls are membranes or semi-transparent mirrors, which keep to/, inside 
but let to/ to escape. We can formally define the membrane operator Ai, which 
separates out a light state (moves it out of the box) and does nothing to a heavy 
state, M.(aj \rrih) + aj |to;}) = (at \m,h) + aj |to/)). Similarly, scattering is described 
by V-operator, V(aj |to;)) = (a' h \mh) + a>'f |to;)), where up and down arrows label the 
amplitudes localized outside and inside the box, respectively. Then the evaporation 
process reads as "separation-scattering-separation", or and 

MVM(aj \m h ) + aj \mi)) = a' h l \m h ) + (aj + af) |m,) . (4) 

The process shown in Fig. [l] can be represented as an operator (A4V) n Xi = 
(MV) . . . (A4V)M., where n is the number of scatterings. This formalism allows 
an easy calculation of the amplitudes, a'' h ^ l = ai™. , as a function of n or time. In 

particular, l&jj 00 '*! 2 = and |aj°°^| 2 = 1, that is, a non-decaying particle can be entirely 
converted into a light mass state and become unconfined. 

Note that scatterings in the box shall not happen too frequently, because the 
mass states have to separate; otherwise they interact coherently as a well-defined 
flavor state. For completeness, we qualitatively consider two more cases. In case (A), 
let's take a box that is completely sealed (no membranes), so both mass states are 
kept inside. Then an equilibrium shall exist and two regimes are possible. First, if 
the wave packets are wide enough (it depends on production/detection, reflection and 
scattering), then they overlap and scatter coherently, hence transition probabilities 
shall be oscillatory functions of space and/or time. Therefore, equilibrium can depend 
on masses and four-momenta of mass states, the box size, L^ OXl and even parameters of 
the scattering potential. Second, when the overlap of wave packets vanishes, scattering 
of mass states occurs independently, so an equilibrium is set by the detailed balance, 
\rrih) <^ | to,;). Equality of the rates of the forward and reverse processes determines 
0^ — the mass and flavor composition of the system. In case (B), we put our original 
"particle-in-a-box" system into a larger sealed box of size Lbox,2 > £box and consider 
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the incoherent regime only (note that light states can come back into the smaller box, 
but heavy states cannot go out of it). The mass conversion is now skewed toward 
\rrih) — > \mi) process, because the rate of the opposite process is reduced by a factor 
of (ibox/ibox,2) 1//2 (since the wave functions are normalized by box volumes). Hence, 
the composition is also skewed towards the light state. In the limit of £box.2 — ► oo, 
one recovers la^ 00 ^ 2 = 0. 

4. Implications 

Possible important implications mentioned in Introduction are estimated here. A more 
detailed analysis is beyond the scope of the paper. 

First, cosmological neutrinos of the CNB. Assuming a nominal non-relativistic 
neutrino cross-section (TU] a v ~ 10~ 60 cm 2 , the mean density of nucleons in halos 
n ~ 10 -3 cm -3 , and the random velocity of v ~ 10 -3 , we estimate the number of 
scatterings to be ~ na v vtn ~ 10~ 47 per neutrino for the life-time of the universe 
(the Hubble time), tjj = 4 x 10 17 s. The effect is very tiny and unobservable (unless 
nonrclativistic neutrinos (self-)interact more strongly than we currently think), though 
it may have some effect for the ultimate fate of the Universe, at tu — > oo. However, 
a v is tremendously enhanced by ~ N 2 Z 2 for coherent scattering, where Z is the 
charge of the atomic nucleus and N ~ nA 3 is the number of nuclei within volume 
~ A 3 and A = 1/mv is the neutrino wavelength |10j . Assuming A ~ 1 cm, for 
Earth (n ~ 10 23 cm~ 3 , nominal Z ~ 25 and distance L ~ 10 4 km) the conversion 
"optical depth" is r c ~ 10~ 4 , that is CNB passing through Earth is slightly affected 
by conversions. 

Some models of dark matter involve a hypothetic sterile neutrino, which is 
presumably more massive than the known ones. The results of the evaporation model 
can be readily generalized to the sterile neutrino scenario. 

In the axion CDM model, axions can scatter off cosmic magnetic 
fields. Using the experimental upper limit on axion-photon coupling g ai ~ 
10 -11 GeV , we estimate the interaction/conversion probability [2 [3] to be ~ 
0.1(g a7 /10~ n GeV _1 )(£?/l juG)(L/10 kpc). For galactic halos with a typical galactic 
field strength of, say, B ~ I — 3 fiG and the halo sizes L ~ 30 — 100 kpc, the proba- 
bility is of order unity per passage. In galaxy clusters with B ~ 0.1 /xG and sizes of 
a couple of Mpc, the rate can be similar. Note that 'one-per-passage' is an efficient 
regime of conversions, so the cosmological effect can be noticeable. If occurs, it can 
lead to substantial evaporation of axion CDM halos on the Hubble time-scale. If the 
effect is not observed, it can further constrain the coupling constant, g aj , or rule out 
axion-CDM. 

The WIMP CDM model deals with the lightest supersymmetric particle 
(ncutralino), higher-mass ones are assumed to decay quickly. However, for high mass 
degeneracy, which we consider now, all decay channels can be kinematically forbidden, 
so more than one particle can be stable. Taking rather crudely Am ~ MeV and 
m ~ TeV, we estimate the velocity of a produced light state to be v', ~ (Am/m) 1 / 2 ~ 
300 km s j which is a cosmologically interesting number. It is smaller or comparable 
to the escape velocity from halos of large galaxies and galaxy clusters, but it is larger 
then v esc for halos of dwarf galaxies, which can be as low as a few tens km s . A halo 
dichotomy is then expected: the number of small halos must be smaller than what 
CDM simulations predict whereas the population of large halos is mainly unaffected. 
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Such an effect can be an alternative explanation to the 'missing satellite problem' and, 
perhaps, it can also affect halo cores [3 [7]. Indeed, mass-conversion is more efficient 
at higher densities (more scatterings). Therefore, higher evaporation rate from high- 
density environments can destroy or smear out small high-concentration halos and 
smooth out the mass distribution in the cores. 

For the effect of conversions to be significant, the rate shall be at least a few 
reactions per Hubble time. Dark matter self-interaction can accommodate the needed 
large cross section, though simple models seem to be disfavored [SJ [7J . The models 
with the velocity-dependent cross-section and the Sommerfeld enhancement [5] , which 
is the most profound at low velocities (that is, in dwarfs, again), shall be tested 
before this scenario is ruled out entirely. Meanwhile, a nice feature of our model 
is that elastic cross-sections for self-interacting CDM and for CDM scattering off 
normal matter, M, can be made as small as desired (to comply with observational 
constraints), while keeping the conversion cross-section large. For instance, for 

scattering h + M ->• h + M and I + M ->• I + M, for 9 ~ tt/4 and V a Vp, 

one readily has Vhh ~ Vu ~ 0, Vih = Vu ~ V a . The disadvantage of this scenario is 
that it required strong degeneracy, which needs an explanation. However, a similar but 
ad hoc model with GeV-scale WIMPs with a keV-scale degeneracy seems to explain 
some direct detection experimental data and spacecraft flyby anomalies [ITJ [T2] . 

To conclude, we have demonstrated that elementary particles with flavor mixing, 
like neutrino, kaons, etc., that are initially trapped in the external gravitational field 
can gradually escape from it in the process involving conversions of mass states. 
In retrospect, the process has a simple interpretation. In ordinary non-relativistic 
quantum physics one typically chooses the zero of energy such that bound states 
have negative energy. In relativistic quantum physics it is more correct to say that 
bound states correspond to quantum states with an energy that is less than the rest 
energy of the particle in question (assuming the potential asymptotes to zero at large 
distances). In a flavour-mixed system, if the energy of the "bound state" is less than 
the rest energy of the initial mass eigenstate the system is prepared in, but greater 
than the rest-energy of any one of the other mass eigenstates the particle can mix 
into, then there will be some probability of escape. For the escape process to proceed 
one needs a flavor-mixed system to experience interactions in the basis that is not 
aligned with the propagation (mass) basis. In this paper we also considered a number 
of possible implications of the conversion effect for the cosmic neutrino background 
and dark matter. We also argue that if dark matter in the universe is a flavor-mixed 
particle, the suggested evaporation process can resolve several outstanding problems 
of the modern cosmology related to the large-scale structure formation of the Universe. 
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